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Abstract 

In this lecture, I explain the gauge-invariant formulation for perturbations of back- 
ground spacetimes with untwisted homologous Einstein fibres, which include lots of 
practically important spacetimes such as static black holes, static black branes and 
rotating black holes in various dimensions. As applications, we discuss the stability of 
static black holes in higher dimensions and flat black branes. 
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§1. Introduction 



Perturbation analysis is a very powerful tool to investigate the dynamical response of a 
system against small disturbances. In particular, in general relativity whose fundamental 
equations are quite hard to solve analytically in general due to their nonlinearity and strong 
couplings, perturbation analysis of exact solutions plays crucial roles in physical and astro- 
physical problems. The most successful example is the perturbative studies of cosmological 
perturbations, which has in particular provided the foundation for the present structure for- 
mation theory and the precise observational cosmology in terms of CMB and gravitational 
waves. 

Another important example is the perturbative studies of black holes. Such an investiga- 
tion was first systematically done for the Schwarzschild black hole by Regge and Wheeler 1 ) in 
1957. In particular, they succeeded in reducing the Einstein equations for odd-parity pertur- 
bations, which is called vector perturbations in the present lecture, to a single master ODE, 
which is called the Regge- Wheeler equation now. The formulation was extended to even- 
parity perturbations (scalar perturbations in this lecture) by Zerilli thirteen years later, 2 ) 
and the master equation called the Zerilli equation was derived for such perturbations. Soon 
later, Teukolsky succeeded in deriving similar Master equations for perturbations of the Kerr 
black hole. 3) 

The original purpose of these formulations appears to have been to study gravitational 
emissions from particles plunging into or orbiting around black holes. However, it was soon 
recognised 4 )' 5 ^ that the formulation can be used to study the stability of black holes, which 
is also a practically important problem in determining the final fate of gravitational collapse. 
Actually, the asymptotically flat neutral and charged black holes were shown to be stable (for 
the proof and its historical background, see the excellent book by Chandrasekhar 6 ) ) . This 
together with the uniqueness theorem for black holes in the asymptotically flat electrovac 
system 7 ) now provide the basis of the current black hole astrophysics. 

These results on four-dimensional black holes are practically sufficient for investigations 
of low energy phenomena. However, taking account of the higher-dimensionality of the 
present candidates for the unified theory, it is likely that higher- dimensional black holes are 
formed in the early universe and in extremely high energy astrophysical phenomena as well 
as in particle accelerators. In fact, motivated by this expectation, lots of work has been done 
on higher-dimensional black holes in various theories, and astonishing discoveries have been 
obtained. In particular, it is now widely recognised that black hole uniqueness does not hold 
in higher dimensions, except for static black holes. 8 )' 9 ) Furthermore, a full list of regular 
black holes has not been obtained even in five dimensions 10 )' 11 ) (Cf. 12), 13)). 
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In this situation, perturbative analysis of exact solutions found so far is expected to be 
quite useful for the study of stability and uniqueness of higher-dimensional black holes. 14 ) In 
particular, it will be a great help if we can reduce the Einstein equations for perturbations 
of higher dimensional black holes to decoupled master equations, as in the four-dimensional 
case. In this lecture, we show that we can really reduce the perturbation equations to 
decoupled mater equations for some classes of black holes and study the stability with the 
help of them. We also point out that such a reduction is not always possible. 

The remaining part is organized as follows. First, in the next section, we briefly overview 
the present status of the black hole stability issue in four and higher dimensions. Then, in 
£j3j we explain the basic aspects of the gauge-invariant formulation for perturbations of a 
general class of background spacetimes that can be written as a warped product of a lower 
dimensional spacetime and an Einstein space. In §H we apply this formulation to static black 
holes and discuss their stability. Next, in §|5j we study the stability of flat black branes and 
point out the non-hermitian nature of the perturbation equations of this system. Section [6] 
is devoted to brief summary and discussion. 

§2. Present Status of the Black Hole Stability Issue 

In this section, we briefly overview the present status of the investigations of the stability 
problem of black holes. It is far from complete. 

2.1. Four Dimensions 

The present status of the stability issue for four-dimensional black holes is summarised 
as follows: 

• Stable 

— Schwarzschild black hole 4 -'' 5 )' 15 )' 16 ) 

— Reissner-Nordstrom black hole 6 ) 

— AdS/dS (charged) black holes 17 )- 18 ) 

— Kerr black hole 19 ) 

— Skyrme black hole (non- unique system) 20 )~ 22 ) 

• Unstable 

— YM black hole (non-unique system) 23 )' 24 ) 

— Kerr-AdS black hole (£f2 h < l,r h < £) 25) 

As is seen from this list, the stability is established for all AF black holes with connected 
horizon in the Einstein-Maxwell system, except for the charged rotating black hole (Kerr- 
Newman black hole). This is because the perturbation equations have not been reduced to 
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decoupled single master equations for this system. 

In the asymptotically adS / dS case, the stability of static black holes have been established 
with the help of mater equations derived by Cardoso and Lentos. 26 -* In contrast, in the 
rotating case, it was conjectured that large Kerr-AdS black holes are stable, while small 
ones are superradiant unstable. 27 )> 28 ) Recently, the conjecture was shown to be true in the 
limit of slow rotation and small horizon. 

2.2. Higher Dimensions 

In contrast to the four-dimensional case, in addition to conventional black holes, there 
exist different kinds of black objects such as black strings, black branes, Kaluza-Klein black 
holes, Kaluza-Klein bubbles and black tubes in higher dimensions. The classification of these 
black objects is far from complete, and rather a little is know about the stability of known 
solutions. For example, concerning the asymptotically flat/dS/adS black holes and black 
branes, the present status of the stability issue is summarised as follows: 

• Stable 

— AF vacuum static (Schwarzschild-Tangherlini) 17 ) 

- AF charged static (D — 5,6— 11) 18 )> 29 ) 

- dS vacuum static (D = 5, 6, 7 - 11), dS charged static (D = 5, 6 - 11)17), is), 29) 

- BPS charged black branes (in type II SUGRA) 30 )' 31 ) 

• Unstable 

— AF/adS static black string and AF black branes (non-BPS) 31 -^ 39 ) 

— Rapidly rotating special Kerr-AdS black holes 40 - 1 

In this list, the stability of static black holes in higher dimensions {D > 4) has been 
proved analytically only for D = 5 in the AF/dS charged case and for D = 5, 6 in the dS 
neutral explained in §H The stability in other dimensions up to D — 11 for these 

black holes was proved numerically. 29 ^ It is expected that the same result holds for D > 11 
as well. 

In contrast, in the asymptotically adS case, stability in D > 4 is not certain even for 
neutral static black holes. This is a delicate problem because instability is expected for 
rotating adS black holes 25 - ),27 ^ 41 ^ 43 ' ) (Cf. Ref. 44)). This instability is understood to arise 
from the combination of superradiance due to a rotating black hole and the time-like nature 
of the adS infinity. Some people conjectured that this superradiance also invokes instabilities 
in doubly spinning black rings 45 ) and Kerr black branes of the form Kerr 4 x M p . 42 ) 

The most impressive result about the stability of black objects in higher dimensions is the 
discovery of the Gregory-Laflamme instability of black strings/branes. 32 - 1 Since then, a large 
amount of work has been done on the classification of black holes and black strings/branes 
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in the iS 1 /torus compactified system and their stability. These researches revealed a rich 
structure of the phase diagram for such systems as well as new instabilities (for review, see 
Refs. 10), 46)). However, no clear understanding has been obtained about the origin and 
fate of these instabilities. It is partly because most of the researches were done by numerical 
methods. In §51 we point out some features that may be obstacles against the analytic 
approach. 

Finally, we have to emphasize that very little is known about the stability of asymptot- 
ically flat solutions with rotating black objects. For example, rapidly rotating Myers-Perry 
solutions were conjectured to suffer from a Gregory-Lanamme type instability, 47 -* but it has 
not been proved by any exact analysis. Black rings are also expected to be unstable because 
of their similarity to black string solutions, but no exact proof has been presented. 

§3. Gauge-invariant Perturbation Theory 

In this section, we explain the basic idea and techniques of the gauge-invariant formu- 
lation of perturbations 48 -*' 49 ) for a class of background spacetimes that includes static black 
hole spacetimes as special case. 

3.1. Background Solution 
3.1.1. Ansatz 

We assume that a background spacetime can be locally written as the warped product 
of a m-dimensional spacetime JV and an n-dimensional Einstein space as 

M n+m ^J/ xX 3 ( Z M) = ^o^ ( 3 .-q 

and has the metric 

ds 2 = g MN dz M dz N = g ab (y)dy a dy b + r(y) 2 da 2 n , (3-2) 

where da 2 = r jijdx l dx^ is an n-dimensional Einstein space J(f satisfying the condition 

R i3 = {n- l)K liy (3-3) 

Note that for n < 3 the Einstein space is automatically a constant curvature space, while 
for n > 3, does not have a constant curvature generically. 

For this type of spacetimes, we can express the covariant derivative Vm, the connec- 
tion coefficients r^ L (z) and the curvature tensor Rmnls(z) in terms of the corresponding 
quantities for Jf m and X n . We denote them D a , m r b a c (y), m R abcd (y) and D u r] k (x), R ijM (x) 
respectively. In particular, the curvature tensor can be expressed as 

D a D b r . . 2 • • 

R a bcd = m R a bcd, R\jb = R l jki = m Rabcd — (Dr) (S l k jji — Sljjk). (3-4) 
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Hence,the non-vanishing components of the Einstein tensor are given by 



Gab 



n 

m G ab - -D a D b r - 
r 



n(n-l)K- (Dr) 



-Dr 
r 



_l mR _ {n-l){n-2)K-{Drf + n-l & 



9- 



(3-5a) 
(3-5b) 



From this and the Einstein equations Gmn + AgMN = k 2 Tmn, it follows that the energy- 
momentum tensor of the background solution should take the form 



T ab = T ab (y) } T ai = 0, Tj = P(y)S i j . 



(3-6) 



3.1.2. Examples 

This class of background spacetimes include quite a large variety of important solutions 
to the Einstein equations in four and higher dimensions. 

1. Robertson- Walker universe: m = 1 and is a constant curvature space. 

ds 2 = -dt 2 + a{t) 2 da 2 n . 

The gauge-invariant formulation was first introduced for perturbations of this back- 
ground by Bardeen 48 ) and applied to realistic cosmological models by the author. 49 ) -51 ) 

2. Braneworld model: m = 2 (and is a constant curvature space). For example, 
the metric of AdS n+2 spacetime can be written 

dr 2 



ds' 



[l-\r l )d1f + r l dQ l n . 



(3-7) 



1 - Xr 2 

The gauge-invariant formulation of this background was first discussed by Muko- 
hyama 52 - 1 and then applied to the braneworld model taking account of the junction 
conditions by the author and collaborators. 53 ) 

Higher-dimensional static Einstein black holes: m = 2 and is a compact 
Einstein space. For example, for the Schwarzschild-Tangherlini black hole, = S n . 
In general, the generalised Birkhoff theorem says 18 ) that the electrovac solutions of the 
form (13-21) with m = 2 to the Einstein equations are exhausted by the Nariai-type 
solutions such that M is the direct product of a two-dimensional constant curvature 
spacetime jV and an Einstein space J?T with r = const and the black hole type solution 
whose metric is given by 

dr 2 



ds 



f(r) 
f(r) = K 



- f(r)dt 2 + r 2 da, 



2M 



-n—l 



+ 



Q 2 



,2n-2 



(3-8) 
(3-9) 
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The gauge-invariant formulation for perturbations was applied to this background to 
discuss the stability of static black holes by the author and collaborators. 17 '*' 18 )' 54 ) This 
application is explained in the next section. 

4. Black branes: m = 2 + k and = Einstein space. In this case, the spacetime factor 
jY is the product of a two-dimensional black hole sector and a fc-dimensional brane 
sector: 

dr 

d s 2 = ^—- f( r )dt 2 + dz-dz + r 2 da 2 . (3-10) 

fin 

One can also generalise this background to introducing a warp factor in front of the 
black hole metric part. The stability of this background for the case in which J?T is an 
Euclidean space is discussed in §51 

5. Higher-dimensional rotating black hole (a special Myers- Perry solution): m = 4 
and JT = S n . 

ds 2 = g rr dr 2 + ggodQ 2 + gudt 2 + 2g t( f > dtd(j) + g^dcf) 2 + r 2 cos 2 9df2 2 , (3T1) 

where all the metric coefficients are functions only of r and 9. The stability of this 
background was studied in Ref. 43). 

6. Axisymmetric spacetime: m is general and n — 1. 

3.2. Perturbations 

3.2.1. Perturbation equations 

In order to describe the spacetime structure and matter configuration (M, g, <P) as a 
perturbation from a fixed background (M,g,<P), we introduce a mapping 

F : background M -> M, (3-12) 

and define perturbation variables on the fixed background spacetime as follows: 

h:=5g = F*g-g, (3-13a) 
<P := 6$ = F*$ - $. (3-13b) 

Then, if the perturbation variables have small amplitudes, the Einstein equations and the 
other equations for matter can be described by linearised equations well. For example, in 
terms of the variable 

= V _ 2 h 9^i ( 3 ' 14 ) 
the linearised Einstein equations can be written as 



2k 2 5T^. (3-15) 
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Gauge 




Background Perturbed 
spacetime spacetime 

Fig. 1. Gauge transformation 
where is the Lichnerowicz operator defined by 

A L Vv := -DVv + + Rued? - 2iW/# a/3 - (3-16) 

3.2.2. Gauge problem 

For a different mapping F', these perturbation variables defined above change their val- 
ues, which has no physical meaning and can be regarded as a kind of gauge freedom. Because 
F and F' are related by a diffeomorphism, the corresponding changes of the variables are 
identical to the transformation of the variables with respect to the transformation / = F'~ 1 F. 
In the framework of linear perturbation theory, we can restrict considerations to infinitesimal 
changes of F. Hence, / is expressed in terms of an infinitesimal transformation £ M as 

6x» = x^f(p))-x^p)=e, (3-17) 
and the gauge transformations are expressed as 

Sh^ = -££9nv = -V^ - V„£ M , (3-18a) 
5cj)=-£ i $. (3-18b) 

From its origin, the perturbation equations including the linearised Einstein equations given 
above are invariant under this gauge transformation. 

To be specific, for our background spacetime, the metric perturbation transforms as 

5h ab = -D a £ b - D b £ a , (3-19a) 
5h ai = -r 2 D a (^) - ACa, (3-19b) 



Shij = -Dfo - D& - 2rD a r£ alij (3- 19c) 
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and the perturbation of the energy-momentum tensor r ; 



8T av transforms as 





—i c D c T ab — T ac D b ^ c — T bc D a ^ c , 
-T ah D£ b -r 2 PD a (r- 2 Z t ), 
-?D a {r 2 P) lij -P{b£ j + b j & 



(3-20a) 
(3-20b) 
(3-20c) 



In order to remove this gauge freedom, one of the following two approaches is adopted 
in general: 

i) Gauge fixing method: this method is direct, but it is rather difficult to find relations 
between perturbation variables in different gauges in general. 

ii) Gauge-invariant method: this method describe the theory only in terms of gauge- 
invariant quantities. Such quantities have non-local expressions in terms of the original 
perturbation variables in general. 

These two approaches are mathematically equivalent, and a gauge- invariant variable can be 
regarded as some perturbation variable in some special gauge in general. Therefore, the 
non-locality of the gauge-invariant variables implies that the relation of two different gauges 
are non-local. 

3.2.3. Tensorial decomposition of perturbations 

In this lecture, we focus on the gauge-invariant approach to perturbations and explain 
that in the class of background spacetimes described above, we can locally construct fun- 
damental gauge invariant variables with help of harmonic expansions. This construction 
becomes more transparent if we decompose the perturbation variables into components of 
specific tensorial types. This decomposition also helps us to divide the coupled set of per- 
turbation equations into decoupled smaller subsets, and in some cases into single master 
equations. 

First of all, note that the basic perturbation variables Hmn and tmn can be classified 
into the following three algebraic types according to their transformation property as tensors 
on the n-dimensional space JtT: 

i) Spatial scalar: h ab , r ab 

ii) Spatial vector: h ai ,rf 

iii) Spatial tensor: ^,rj 

Among these, spatial vectors and tensors can be further decomposed into more basic 
quantities. First, we decompose a vector field Vi on Jtf into a scalar field and a transverse 
vector v!P as 



(3-21) 
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Then, from the relation 

Av {s) = Div\ (3-22) 

the component fields and vf^ can be uniquely determined from Vi up to the ineffective 
freedom in v ^ to add a constant, provided that this Poisson equation has a unique solution 
on up to the same freedom. For example, when is compact and closed, this condition 
is satisfied. 

Next, we decompose a symmetric tensor field of rank 2 on as 

Uj = -tgtj + DiDjS - -As gij + Drfj + DjU + tf^; (3-23a) 
/).f = 0, 4 tt}1 = 0, Dit { - t)[ = 0. (3-23b) 

Here, t is uniquely determined as t = t\. Further, from the relations derived from this 
definition, 

A(A + nK)s = — ^— \bibd ij - -At] , (3-24a) 
n — 1 \ n J 

[A + (n - l)K]t l = (5) - b i A~ 1 D j )(b m t jm - n^DH), (3-24b) 

s and ti, hence tf*\ can be uniquely determined from up to the addition of ineffective 
zero modes, provided that these Poisson equations have solutions unique up to the same 
ineffective freedom. 

After these decompositions of vectors and tensors to basic components, we can classify 
these components into the following three types: 

i) Scalar type: v i = £)V s) , = Hg^ + b { bjS - ± Asgij. 

ii) Vector type: V{ = vf\ Uj = D{tj + DjU. 

iii) Tensor type: v l = 0, t) = tf Y . 

We call these types reduced tensorial types. In the linearised Einstein equations, through 
the covariant differentiation and tensor-algebraic operations, quantities of different algebraic 
tensorial types can appear in each equation. However, in the case in which Jtf is a constant 
curvature space, perturbation variables belonging to different reduced tensorial types do not 
couple in the linearised Einstein equations if we decompose these perturbation equations into 
reduced tensorial types as well, because there exists no quantity of the vector or the tensor 
type in the background except for the metric tensor. The same result holds even in the case 
in which is an Einstein space with non-constant curvature, because the only non-trivial 
background tensor other than the metric is the Weyl tensor that can only transform a 2nd 
rank tensor to a 2nd rank tensor. 
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Here, note that gauge transformations can be also decomposed into reduced tensorial 
types, and the gauge transformation of each type affects only the decomposed perturba- 
tion variables of the same reduced tensorial type. Hence, gauge-invariant variables can be 
constructed in each reduced tensorial types independently. 

3.3. Tensor Perturbation 

Let us start from the tensor- type perturbation, for which the argument is simplest. 
3.3.1. Tensor Harmonics 

We utilise tensor harmonics to expand tensor-type perturbations. They are defined as 
the basis for 2nd-rank symmetric tensor fields satisfying the following eigen-value problem: 

(A L - A L )T 4J = 0; T] = 0, D,T i= 0, (3-25) 

where A^ is the Lichnerowicz operator on defined by 

A L hij := —D ■ Dhij - 2R ikjl h kl + 2(n - VjKhy. (3-26) 

When Jtf is a constant curvature space, this operator is related to the Laplace-Beltrami 
operator by 

A L = -A + 2nK, (3-27) 

and, satisfies 

(A + k 2 )Tij = 0; k 2 = X L - 2nK. (3-28) 

We use k 2 in the meaning of — 2nK from now on when JtT is an Einstein space with 
non-constant sectional curvature. 

The harmonic tensor has the following basic properties: 

1. Identities: Let T^- be a symmetric tensor of rank 2 satisfying 

T\ = 0, D j Tij = 0. 

Then, the following identities hold: 

2D {i T j]k D [i T j]k = 2D i (T jk D [i T j]k ) + T jk [-AT jk + RjT lk + Ri jk iT a ] , 
2D {i T j)k D (i T j)k = 2D i {T jk D (i T j)k ) + T jk [-AT jk - R\T lk - Ri jk iT il ] . 

On the constant curvature space with sectional curvature K, these identities read 

2D [i T j]k D [i T j]k = 2D i (T jk D [i T j]k ) + T jk (-A + nK)T jk , 
2D {i T j)k D (i T j)k = 2D\T jk D {i T j)k ) + T jk (-A - nK)T jk . 
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2. Spectrum: When is a compact and closed space with constant sectional curvature 
K, these identities lead to the following condition on the spectrum of k 2 : 

k 2 >n\K\. (3-29) 

In contrast, when JfT is not a constant curvature space, no general lower bound on the 
spectrum k 2 is known. 

3. When is a two-dimensional surface with a constant curvature K, a symmetric 2nd 
rank harmonic tensor that is regular everywhere can exist only forfT < 0: for T 2 (K = 
0), the corresponding harmonic tensor T^- becomes a constant tensor in the coordinate 
system such that the metric is written ds 2 = dx 2 + dy 2 (k 2 = 0); for H 2 /T(K = —1), a 
harmonic tensor corresponds to an infinitesimal deformation of the moduli parameters. 

4. For Jtf = S n , the spectrum of k 2 is given by 

A; 2 = /(/ + n-l)-2; I = 2, 3, • • • , (3-30) 

3.3.2. Perturbation equations 

The metric and energy-momentum perturbations can be expanded in terms of the tensor 
harmonics as 

h ab = 0, h ai = 0, h ij = 2r 2 H T T ij , (3-31) 

r ab = 0, 7f = 0, t) = t t T). (3-32) 

Since the coordinate transformations contains no tensor- type component, Ht and tt are 
gauge invariant by themselves: 

£ M = Sz M = 0; (3-33a) 
5H T = 0, 5t t = 0. (3-33b) 

Only the (i, j)-component of the Einstein equations has the tensor-type component: 

77 k 2 4- IK 
- DH T - -Dr ■ DH T + —~ 2 H T = R 2 t t . (3-34) 

Here, □ = D a D a is the D'Alembertian in the m-dimensional spacetime jY . Thus, the 
Einstein equations for tensor-type perturbations can be always reduced to the single master 
equation on our background spacetime. 

3.4. Vector Perturbation 

3.4.1. Vector harmonics 

We expand transverse vector fields in terms of the complete set of harmonic vectors 
defined by the eigenvalue problem 

(A + k 2 )Yi = 0; AV* = 0. (3-35) 
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Tensor fields of the vector-type can be expanded in terms of the harmonic tensors derived 
from these vector harmonics as 

V« = ~(AV i + 4-V i ). (3-36) 

They satisfy 

A + k 2 - (n + 1)k\ Y i:j = 0, (3-37a) 

Vj = 0, DM = {n 2 ~ 1)K V„ (3-37b) 

Here, there is one subtle point; vanishes when Vj is a Killing vector. For this mode, 
from the above relations, we have 

k 2 = {n - l)K. (3-38) 

We will see below that the converse holds when J{f is compact and closed. We call these 
modes exceptional modes. 

Now, we list up some basic properties of the vector harmonic relevant to the subsequent 
discussions. 

1. Spectrum: In an n-dimensional Einstein space satisfying 

Rij = {n- l)K 9ij , (3-39) 

we have 

ID^^V^ = 2Di(yjD^V j] ) + Vj [-A + (n - 1)K] V j , (3-40a) 
2D (i ^ ) £)( i y J ' ) = 2D i {V j D^ i V i) ) + Vj [-A - (n - 1)K] V j . (340b) 

When is compact and closed, from the integration of these over JT, we obtain the 
following general restriction on the spectrum of k 2 : 

k 2 > (n — 1)\K\. (3-41) 

Here, when the equality holds, the corresponding harmonic vector becomes a Killing 
vector for K > and a harmonic 1-form for K < 0, respectively. 

2. ForJT" = S n , we have 

k 2 = £(£ + n - 1) - 1, (£=1,2,---)- (3-42) 

Here, the harmonic vector field V, becomes a Killing vector for I — 1 and is exceptional. 

3. For K — 0, the exceptional mode exists only when is isometric to T p x ^ n ~p ? where 
^ n_p is a Ricci flat space with no Killing vector. 
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3.4.2. Perturbation equations 

Vector perturbations of the metric and the energy-momentum tensor can be expanded 
in terms of the vector harmonics as 



h ab = 0, h ai = rf a Y i , = 2r 2 H T Y lj , (343a) 
r ab = 0, 7? = rr.V,, rj = r T N). (343b) 

For the vector-type gauge transformation 

i a = 0, 6 = rLW t (3-44) 

the perturbation variables transform as 

6 f a = -rD a (jj , 5H T = *L, 5r a = 0, 5r T = 0. (3-45) 

Hence, we adopt the following combinations as the fundamental gauge-invariant variables 
for the vector perturbation: 

T 

generic modes: T a , r T , F a = f a + j; D aH T (3-46) 
exceptional modes:r a , = rD a ~ r 'Db (J^j (3-47) 

Note that for exceptional modes, F a = f a because H T is not defined. 

The reduced vector part of the Einstein equations come from the components correspond- 
ing to G\ and G*. In terms of the gauge-invariant variables defined above, these equations 
can be written as follows. 

• Generic modes: 

1 D b (r n+1 F {1) ^ F ( " ' A 



r n+l \ ab / r 



\ k 2 — (n — 1)K 2 
) 72 — F - = - 2Rt -> 



(348a) 

^-D a (r n - 1 F a ) = -R 2 r T . (348b) 

Exceptional modes: k 2 = (n — \)K > 0. For these modes, the second of the above 
equations coming from does not exist. 

^-D b (r^Fff) = -2R 2 r a . (3-49) 
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3.5. Scalar Perturbation 
3.5.1. Scalar harmonics 

Scalar functions on can be expanded in terms of the harmonic functions defined by 

(A + A; 2 )§ = 0. (3-50) 

Correspondingly, scalar-type vector and tensor fields can be expanded in terms of harmonic 
vectors §j and harmonic tensors §y define by 

Si = ~Di$, (3-51a) 
k 

S ij = -^A4- S +-7«S. (3-51b) 
These harmonic tensors satisfy the following relations: 

D& = fcS, (3-52a) 
[A + k 2 - (n - l)K]§i = 0, (3-52b) 



n — 1 k 2 — nK . 
n k 

[A + k 2 - 2nK]§ ij = 0. (3-52d) 



0, £)$ = , Si, (3-52c) 



Note that as in the case of vector harmonics, there are some exceptional modes: 

i) k = 0: Sj = 0, % = 0. 

ii) k 2 =nK (K > 0): §^ = 0. 

For scalar harmonics, k 2 = is obviously always the allowed lowest eigenvalue. Therefore, 
the information on the second eigenvalue is important. In general, it is difficult to find such 
information. However, when J%f n is a compact Einstein space with K > 0, we can obtain a 
useful constraint as follows. Let us define Qij by 

Qij : DiDjY - - gij AY. 

n 



Then, we have the identity 

QijQ ij = D i (D i Y DiDjY - YD { AY - R i:j D j Y) + Y [A(A + (n - l)K)\ Y - —(AY)' 
For Y — S, integrating this identity, we obtain the constraint on the second eigenvalue 

k 2 > nK. (3-53) 
For = S n , the equality holds because the full spectrum is given by 

k 2 = £(£ + n -l), (£ = 0,1,2, ■■■). (3-54) 
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3.5.2. Perturbation equations 

The scalar perturbation of the metric and the energy-momentum tensor can be expanded 

as 



h ab = /ofcS, h ai = rf a 8>i, hij = 2r 2 (H L -f i:j S + H T E>ij), (3-55a) 

Tab = r ab §, if = rr a Si, rj = 5P5fi + r T S). (3-55b) 

For the scalar-type gauge transformation 

U = T a §, & = rLS i: (3-56) 

these harmonic expansion coefficients for generic modes k 2 (k 2 — nK) > of a scalar- type 
perturbation transform as 

6 U = -D a T b - D b T a , (3-57a) 

6 f a = -rD a (^j + k -T a , (3-57b) 
k D a r 

SH L = L T a , (3-57c) 

nr r 

5H T = -L, (3-57d) 

r 

5rab = —T c D c T ab - f ac D b T c - T bc D a T c , (3-57e) 

5r a = -{TabT h -PT a ), (3-57f) 

r 

5(5P) = -T a D a P, (3-57g) 

5r T = 0. (3-57h) 

From these we obtain 

5X a = T a ; X a = T - (/ a + £a,#t) • (3-58) 

Hence, the fundamental gauge invariants can be given by tt and the following combinations: 

F = H L + -H T + -D a rX a , (3-59a) 
n r 

Fab = fab + D a X b + D b X a , (3-59b) 

^ = r ab + f b c D a X c + f c a D b X c + X c D c f ab , (3-59c) 

= r a - ^(f a % - PX a ), (3-59d) 

r L = <5P + X a D a P. (3-59e) 
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The scalar part of the Einstein equations comes from G a b, G a i and Gj. First, from SG a b, 
we obtain 

-UF ab + D a D c F b c + D b D c F c a + n ^(-D c F ab + D a F cb + D b F ca ) 

( k 2 - -\ 

+ m R c a F cb + m R c b F ca — 2 m R acbd F c + — — R + 2A\ F ab — D a D b F^ 



■In | /A,/;,/'' + -D a rD b F + -D b rD a F 

r r 



D c D d F cd + —D c rD d F cd + (—D c D d r + n ^ - l) D c rD d r 

/y \ 7" 7" 



mrycd 



)F cd - 



2nDF- 



2n 
r 

2n(n + 1 



h 2 — nK 

Dr-DF + 2(n-l)- r— — F 
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-□F c c - -L>r • DF C C + — F c c 



9ab — 2n 2 E ab . 



Second, from we obtain 



k 
r 



-D b (r n - 2 F b a ) + rD a [ -F b ) + 2(n - l)D a F 



— 2/t S n 



Finally, from the trace- free part of SG l j, we obtain 

k 2 



and from the trace 5G], 



--D a D b F ab - 
2 r 



2r 2 

n — 1 



[2( n -2)F + F a a ]=«: 2 r T , 



+ 



2 



2r 2 



D a rD b F ab 

D a rD b r - (n - 1)- 



(n _l) (n _ 2 ) ^D»D*r 



-Dr ■ DF C C - ^-^%F C C + (n - l)DF 
2n r 2 



+ ^21 D r • DF - '"- 1 ""- 2 » t2 -" A F = #E 



(3-60) 



(3-61) 



(3-62) 



L- (3-63) 

r n 

Note that for the exceptional mode with k 2 = nK > 0, the third equation does not exist, 
and for the mode with k 2 = 0, the second and the third equations do not exist. For these 
exceptional modes, the other equations hold without change, but the variables introduced 
in the above are not gauge invariant. 

Although the energy- momentum conservation equation VatT^ = can be derived from 
the Einstein equations, it is often useful to know its explicit form. For scalar-type perturba- 
tions, they are given by the following two sets of equations: 



n 



kr 
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k 

2r' 



+ ^-(T ab F ab 



PF„ a ) = 0, 



+T b D b F - PD a F + i (T*D b F° - T bc D a F bc ) 
§4. Stability of Static Black Holes 



0. 



(3-64a) 



(3-64b) 



We study the stability of static black holes utilising the gauge-invariant formulation 
for perturbations explained in the previous section. We consider the static Einstein black 
hole which corresponds to the case with m = 2 of the general background considered in the 
previous section and has the metric (I3-9P . The key point is the fact that gauge-invariant per- 
turbation equations can be reduced to decoupled single master equations of the Schrodinger 
type for any type of perturbations in this background. 

4.1. Tensor Perturbations 

The gauge-invariant equation for tensor perturbations is already given by a single equa- 
tion for each mode, Assuming that the source term vanishes, it reads 



d t H 2 +fd r (fd r H T ) 



k 2 + 2K 



fH T = 0. 



(4-1) 



Here, note that even if there exist electromagnetic fields, tj- vanishes because the electro- 
magnetic field is vector-like and does not produce a tensor-type quantity in the linear order 
at least. 

With the help of the Fourier transformation with respect to t, i.e., assuming Hq^ oc c %u) , 
this equation can be put into the Schrodinger-type eigenvalue problem; 

H T = r~ nl2 $(r)e- iuJ \ 



-fd r (fd r 0) + V t <P 



where 



f 



f 



k 2 + 2K + 



+ 



k 2 + 



n — 2n + 



nrf n{n — 2)f 
"1 

n(n + 2) 



(4-2) 
(4-3) 

(4-4) 



-K 



Xr 2 + 



n 2 M n(3n - 2)Q 5 



I | 2r n ~ 1 4 r 2n-2 

If Vt is non-negative, we can directly conclude the stability. However, it is not so easy to 
see whether V t is non-negative or not outside the horizon. This technical difficulty is easily 
resolved by considering the energy integral 



E :-- 



dr 



-Ad t H T ) 2 + f(d r H T ) 2 + k -^E 2 T 



(4-5) 
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From the equation for Ht, we find that 

d t E = 2 [fd t H T d r H T ] r ~ = 0. (4-6) 

Hence, in the case is a constant curvature space, the condition on the spectrum k 2 > n\K\ 
guarantees the positivity of all terms in E, and as a consequence the stability of the system. 

4.2. Vector Perturbations 
4.2.1. Master equation 

For vector perturbations, the energy-momentum conservation law is written 

777 

D a (r n+1 r a ) + -^r n r T = 0. (4-7) 

For m v = k 2 — (n — 1)K ^ 0, with the help of this equation, the second of the perturbation 
equations, (13-48bl) . can be written 

2k 2 

D a (r n - l F a ) = D a {r n+l r a ). (4-8) 

ui v 

In the case of m = 2, from this it follows that F a can be written in terms of a variable Q as 



2k 2 

r n~lpa = e ab DbQ + r ^ T \ (4.9) 

m v 

Further, the first of the perturbation equations, ( l3-48ah . is equivalent to 

D a (r n+1 F^) - m v r n ~ l e ab F h = -2K 2 r n+1 e ab r b , (4-10) 
where e a b is the two-dimensional Levi-Civita tensor for g abl and 

F« = e ab rD a = e ab rD a H±\ . (4-11) 



Inserting the expression for F a in terms of Q into (14-101) . we obtain the master equation 

L D a Q \ -^Q = - 2 J^r n e ab D a {rr h ). (4-12) 
r n J r z m v 

Next, for m v = 0, the perturbation variables Ht and tt do not exist. The matter 
variable r a is still gauge- invariant, but concerning the metric variables, only the combination 
F^ defined in terms of f a in (14- lip is gauge invariant. In this case, the Einstein equations 
are reduced to the single equation (14-101) . and the energy-momentum conservation law is 
given by (14-71) without the tt term. Hence, r a can be expressed in terms of a function 
as 

r n+1 T a = e ab D b r^. (4-13) 
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Inserting this expression into (14401) with e cd D c (F d /r) replaced by F (1) /V, we obtain 

D a (r n+1 F^) = -2n 2 D a r {l \ (4-14) 

Taking account of the freedom of adding a constant in the definition of the general 
solution can be written 



2/tM 1 ) 



(4-15) 



Hence, there exists no dynamical freedom in these special modes. In particular, in the 
source-free case in which is a constant and K — 1, this solution corresponds to adding 
a small rotation to the background static black hole solution. 

4.2.2. Neutral black holes 

For a neutral static Einstein black hole, the master equation for a generic mode can be 
put into the canonical form as 

u 2 $=-fd r (fd r $)+V v $ 



(4-16) 
(4-17) 



where 



Vr 



f 



f 



nrf n(n + 2)f 
m v 1 ; 



, n(n + 2)K n(n - 2) 2 3n 2 M 



k z + 



(4-18) 



4 4 2r n - x 

This equation is identical to the Regge- Wheeler equation for n = 2, K = 1 and A = 0. 
In this case, we can put V v into an obviously non-negative form as 



v v = ^ K + 3/) 



(4-19) 



proving the stability of the black hole against vector perturbations (or axial or odd pertur- 
bations). 

In higher dimensions, the potential V v is not positive definite anymore and we can not 
use this type of argument. However, we can still prove the stability with the help of the 
conserved energy integral as in the case of tensor perturbations. In the present case, if we 
define E as 



E :-- 



we have 



E 



-{d t Q) 2 + f{d r Q) 



7 

— d t f2d r f2 



2 + ^ Q 2 



0. 



(4-20) 



(4-21) 



Further, all terms of E is non-negative because m v > 0. Hence, the stability can be con- 
cluded. 
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1=0.05 
K=1, 1=2 





Fig. 2. V v for K = 1,X 
0,1 = 2. 



Fig. 3. V v for K = 1,A > Fig. 4. V v for if = 1,A < 
0, Z = 2. 0, Z = 2. 



4.2.3. Charged black hole 

The formulation for neutral static black holes can be extended to charged static black 
holes. The final master equations consist of two equations: the extension of the equation for 
gravitational perturbations with an electromagnetic source and the equation coming from 
the Maxwell equations: 



m v + 2(n - l)K 



D ( r n-2 D a 



si 



(4-22a) 

'' [m v f2 + 2K 2 q#f] , (4-22b) 



„'2n 



where stf is the gauge invariant representing a vector perturbation of the vector potential of 
the electromagnetic field defined by 



5A a = 0,5A i = ^Y i , 



(4-23) 



and q is the black hole charge related to the charge parameter Q in the background metric 
by 



Q 2 



K 2 q 2 



(4-24) 



n(n — 1) 

By taking appropriate combinations, these equations can be transformed to the two 
decoupled equations 

- = (-C + V±)$±, (4-25) 
where the effective potentials are given by 
/ 



n(n + 2)K n(n - 2) , 2 n(5n - 2)Q 2 /z± 
m v H : : Ar H : — h 



n-1 



(4-26) 
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with 

n 2 4- 9 

l*± = ^— M ± A; (4-27) 

A 2 = (n 2 - l) 2 M 2 + 2n(n - I)m^Q 2 . (4-28) 

4.2.4. S- deformation 

The effective potentials V± are not positive definite as in the neutral case. In the present 
case, we prove that the system is still stable not by the energy integral method, but rather 
by a different method, which we call the S- deformation. 17 ^ 

We first explain the basic idea by the eigen value equation 

.2, 



■<p = [-D 2 + V(r)) $, (4-29) 



where D = d rt . If there exists an unstable mode with uo 2 < 0, we can show $ falls off 
sufficiently rapidly at horizon and at infinity if V is non-negative at horizon and at infinity. 
Hence, we obtain the integral identity, 

u 2 / |<?| 2 ^ = / [|D*| 2 + V(r)\<P\ 2 ] ^. (4-30) 

Jr h J Jr h J 

If V(r) is non- negative definite, this leads to contradiction and hence proves the stability 
because the right-hand side is non- negative. In contrast, in the case in which the sign of V 
is not definite, we cannot say anything about stability from this equation. 

In order to treat such a case, let us replace D by D = D—S. Then, by partial integrations, 
we obtain the modified integral identity with D and V replaced by D and V given by 

V = V + f^-S 2 . (4-31) 
dr 

Hence, if we can find S such that the modified effective potential V is non-negative, we can 
establish the stability of the system even when the original potential is not non-negative 
definite. 

For example, by the S-transformation with 

S = ^t, (4-32) 
the effective potentials V± above can be modified into 



V ± =V ± + ff-S 2 
dr 



3n 



2 



m v + —[—M + ^ ± 



(4-33) 



Here, V + is obviously positive definite. We can also show that V_ is also positive definite. 
Hence, a charged static Einstein black hole is stable for vector perturbations. 
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4.3. Scalar Perturbations 

4.3.1. Master equation 

For a static Einstein black hole background, assuming that F a b, F oc e~* w *, we can reduce 
the whole linearised Einstein equations into a single master equation, as in the case of vector 
perturbations: 54 ^ 

u 2 <P = -fd r (fd r <P)+V s $, (4-34) 
where the master variable <P is defined as 

<P= (2F+-^); (4-35) 

H = m+ n{n r + J M , m = k?-nK, (4-36) 
and the effective potential V s is given by 

y.(r) - (4-37) 
U(r) = - [n 3 (n + 2)(n + 1) V - 12n 2 (n + l)(n - 2)mi 
+4(n - 2)(n - 4)m 2 ] Ar 2 + n 4 (n + 1) V 

+n(n + 1) [4(2n 2 - 3n + 4)m + n(n - 2)(n - 4)(ra + x 2 

-12ra [(n - 4)m + n(n + l)(n - 2)X] tra 

+16m 3 + 4i^n(?i + 2)m 2 , (4-38) 

with x = 2M/r n ~ 1 . 

4.3.2. Neutral black holes 

The above master equation is identical to the Zerilli equation for the four-dimensional 
Schwarzschild black hole (n — 2, K — 1 and A = 0). In this case, from 

/ / 2/ «n 6m2M 36mM 2 72M 3 \ ft , AnnS 
V s = ^[m 2 (m + 2) + — + — — + — — >0, (4-39) 



where m = (/ — 1)(/ + 2) (7 = 2, 3, • • • ), we can easily prove the stability of the black hole. 
In higher dimensions, however, the effective potential V s is not positive definite. Hence, an 
instability may arise. 

Nevertheless, in the case of K — 1 and A = 0, i.e. for the Schwarzschild- Tangherlini 
black hole, we can prove the stability by applying the S-deformation to the energy integral. 
First, from the above master equation, we obtain 

f r °° rlr 

E := — [(d t <P) 2 + (D<P) 2 + V S <P 2 ] , (4-40) 

Jr J 

E = [2fd t $d r <P] r - = 0, (4-41) 
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Fig. 5. V s for K = 1,A = Fig. 6. V s for K = 1, A > Fig. 7. V s for if = 1, A < 
0,Z = 2. 0,Z = 2. 0,Z = 2. 



where .D = fd r . Next, we replace D to 

D = fd r + S. (4-42) 

Then, by partial integration we obtain 



E 

where 

For example, for 



/ - (d t $) 2 + (D$) 2 + V S <P 2 , (4-43) 

Jro J 

V S = V S + f^- - S 2 . (4-44) 



we obtain 



3=—!-, h ^r^ 1 - 1 {(l - + n) + n(n + l)x/2} . (4-45) 

V = (4-46) 

s 4r 2 {(I- + n) + n(n + l)x/2} ' 1 ' 

where 

Q(r) = /x[/n(n + l)x + 2(1 - \){n 2 + n(3Z - 2) + (I - l) 2 }] . (4-47) 

Clearly V s > 0. 

4.3.3. Charged black holes 

For charged black holes, we can also reduce the perturbation equations to decoupled single 
master equations. First, we generalise the master variable <P for the metric perturbation given 
in (14-351) by replacing H by 

H = m + n(n + l)M_rW 

r n—l r 2n—2 v ' 
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Next, we introduce the gauge-invariant variable srf in terms of which the scalar perturbation 
of the electromagnetic field is expressed as 

63? ab + D c {E X c )e ab S, = <fe a6 §, (4-49a) 
S&ai - kE e ab X% = re ab £%, (449b) 
= 0, (449c) 



with 



@a — —^—[D a £rf , r n S = -k 2 ^+UF c c -2nF). (4-50) 



Then, the Einstein and Maxwell equations for scalar perturbations of a charge Einstein 
black hole can be reduced to the following two coupled equations: 18 ^ 



where Vg is the extension to the charged case of the corresponding potential in the neutral 

case, 

_ f(r)U s (r) _ 



U s = [-n 3 (n + 2)(n + l) 2 x 2 + 4n 2 (n + l){n(n 2 + 6n - 4)z + 3(n - 2)m}x 
-12n 5 (3n - 2)z 2 - 8n 2 (lln 2 - 26n + 12)mz - 4(n - 2)(n - 4)m 2 ] Ar 2 
+n 4 (n + 1)V + n(n + 1) {-3n 2 (5n 2 - 5n + 2)z + 4(2n 2 - 3n + 4)m 
+n(n - 2)(n - 4)(n + 1)^} x 2 + 4n [n 2 (4n 3 + 5n 2 - lOn + 4)^ 2 
- {n(34 - 43n + 21n 2 )m + n 2 (n + l)(n 2 - lOn + 12)K} z 
-3(n - 4)m 2 - 3n(n + l)(n - 2) Km] x - 4n 5 (3n - 2)z 3 
+12n 2 {2(-6n + 3n 2 + 4)m + n 2 (3n - 4)(n - 2)X} z 2 
+ {4(13n - 4)(n - 2)m 2 + 8n 2 (lln 2 - 18n + 4)Km] z 

+I6m 3 + 4n(n + 2)Km 2 , (4-53) 

with z = g 2 /r 2 (™- 1 ), and P S1 and P z are the functions of r given by 

P 51 = [-An 4 z + 2n 2 (n + l)x - An{n - 2)m] Ar 2 

+ [2n 2 (n - l)x + An(n - 2)m + 4n 3 (n - 2)K) z 
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-n 2 (n 2 - l)x 2 + {-4n(n - 2)m + 2n 2 (n + 1)K) x 

+4m 2 + 4n 2 mK, (4-54) 
P z = [- n 2 {n + l)x + 2n 2 (3n - 2)z + 2(n - 2)m] Ar 2 

+n(n + l)x 2 + {n 2 (3n - 7)z + (4n - 2)m + n(n + l)(n - 2)fT} x 

-2n 2 (n - 2)^ 2 - {(6n - 4)m + 2n 2 (3n - 4)fT} z - 2nmir. (4-55) 

As in the case of vector perturbations, we can find linear combinations of stf and <£>, in 
terms of which these equations are transformed to the decoupled equations 

7* ± = -(/*±)' + ^* ± ; v ± = w ± > < 4 «) 

Here, 

#+ = 1 y — — -5x, H^ = m + -i— — -(1 + mS)x, (4-57) 

and 5 is a non-negative constant determined from Q by 

Q 2 = (n + l) 2 M 2 5(l + m5). (4-58) 

The effective potentials U± can be expressed in terms of x, Ar 2 , m and 5 as follows: 

C/ + = [-4n 3 (n + 2)(n + 1) 2 5V - 48n 2 (n + l)(n - 2)&r 

-16(n - 2)(n - 4)] Ar 2 - 5 3 n 3 (3n - 2)(n + 1) 4 (1 + m5):r 4 
+4<5 2 n 2 (n + l) 2 {(n + l)(3n - 2)m5 + 4n 2 + n - 2} 

+4<5(n + 1) {(n - 2)(ra - 4)(ra + l)(m + n 2 K)5 - 7n 3 + 7n 2 - 14ra + 8} x 2 

+ {I6(n + 1) (-4m + 3n 2 (n - 2)K) 5 - 16(3n - 2)(n - 2)} x 

+64m + 16ra(n + 2) if, (4-59) 



[/_ = [-4n 3 (n + 2)(n + 1) 2 (1 + mSfx 2 + 48n 2 (n + l)(n - 2)m(l + m5)x 
-16(n - 2)(n - 4)m 2 ] y - n 3 (3n - 2)(n + 1) 4 5(1 + rrwJ) V 
-4n 2 (n + 1) 2 (1 + mS) 2 {(n + l)(3ra - 2)m5 - n 2 } x 3 
+4(n + 1)(1 + mS) {m(n - 2)(n - 4)(n + l)(m + n 2 fT)5 
+4n(2n 2 - 3n + 4)m + n 2 (n - 2)(ra - 4)(n + l)if} a; 2 
-16m {(n + l)m (-4m + 3n 2 {n - 2)X) 5 
+3n(n - 4)m + 3n 2 (n + l)(n - 2)7^} x 

+64m 3 + 16n(n + 2)m 2 K (4-60) 
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Table I. stability of generalised static black holes. 





Tensor 
VQ 


Vector 
VQ 


Scalar 


O — 




K = 1 


A = 


OK 


OK 


OK 


D = 4,5 OK 
£) > 6 ? 


A > 


OK 


OK 


D < 6 OK 
D > 7 ? 


D = 4,5 OK 
D > 6 ? 


A < 


OK 


OK 


D = 4 OK 

D > 5 ? 


D = 40K 

D > 5 ? 


K = 


A < 


OK 


OK 


D = 40K 

D > 5 ? 


D = 40K 

D > 5 ? 


K = -l 


A < 


OK 


OK 


D = 4 0K 
£> > 5 ? 


D = 40K 

D > 5 ? 



By applying the S'-deformation to V + with 

„ / dh 



h + dr 

we obtain 



h+ = r n > 2 ~ 1 H + , (4-61) 



; 2 P 

Vs+ = 2^h; [(n " 2)(n + + 2] • (4 ' 62) 

Since this is positive definite, the electromagnetic mode <P + is always stable for any values 
of K, M, Q and A, provided that the spacetime contains a regular black hole, although V + 
has a negative region near the horizon when A < and Q 2 /M 2 is small. 

Using a similar transformation, we can also prove the stability of the gravitational mode 
<P_ for some special cases. For example, the ^-deformation of V_ with 

s= ±dh_ 
h- dr 

leads to 



h_ = r nt2 - l H_ (4-63) 



7 2 e 

V_ = —t— [2m - (n + l)(n - 2)(1 + mS)x\ . (4-64) 

For n — 2, this is positive definite for m > 0. When X = 1, A > and n = 3 or when 
A > 0, Q = and the horizon is S 4 , from m > n + 2 (/ > 2) and the behaviour of the horizon 
value of x (see Ref. 18 ) for details), we can show that V S - > 0. Hence, in these special cases, 
the black hole is stable with respect to any type of perturbation. 

However, for the other cases, Vs- is not positive definite for generic values of the param- 
eters. The .S-deformation used to prove the stability of neutral black holes is not effective 
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either. Recently, Konoplya and Zhidenko studied the stability of this system for n > 2 
numerically. They found that if A > 0, the system is stable for n < 9, i.e., D > 11. 29 ) 

4.4. Summary of the Stability Analysis 

The results of the stability analysis in this section can be summarised in Table [H In this 
table, D represents the spacetime dimension, n+2. The results for tensor perturbations apply 
only for maximally symmetric black holes, while those for vector and scalar perturbations 
are valid for black holes with generic Einstein horizons, except in the case with K — 1, Q — 
0, A > and D = 6. 

Note that this is a summary of the analytic study. As we mentioned above, the stability 
of AF/dS black hole is shown for D < 12 numerically. Stability is shown also for topological 
adS black holes with non-positive mass. 55 ) 

§5. Flat black brane 

Static flat black brane solutions are perturbatively unstable in contrast to asymptotically 
simple static black holes discussed in the previous section. This was first shown by Gregory 
and Laflamme for the s-mode perturbation, i.e. perturbations that is spherically symmetric 
in the directions perpendicular to the brane. 32 )' 33 - 1 Later on, it was shown that the system 
has no other unstable modes numerically. 37 )> 38 ) These analyses however assumed that the 
frequency of an unstable mode, if it exists, is pure imaginary. In the static system this 
assumption may appear to be natural, but it is not the case in reality. In this section, we 
explain this point explicitly by applying the the gauge-invariant formulation in the previous 
section to this system. 

5.1. Strategy 

Let us rewrite the (m + n + 2)-dimensional flat black brane solution 

ds 2 = -f{r)alt 2 + /(r)-Mr 2 + r 2 da 2 n + dx 2 , (5-1) 

which is the product of (n + 2)-dimensional static black hole solution and the m-dimensional 
Euclidean space, as 

ds 2 = g ab (y)dy a dy b + r 2 da 2 n (5-2) 
with the (m + 2)-dimensional metric 

ds 2 rn+2 = g ab (y)dy a dy b = -f(r)dt 2 + f{r)- l dr 2 + dx 2 . (5-3) 

Then, we can classify metric perturbations into tensor, vector and scalar types with respect 
to the n-dimensional constant curvature space Jif n with the metric da 2 = r )ij{z)dz l dz^ and 
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apply the gauge-invariant formulation developed in the previous chapter to them. Further, 
since the background spacetime is homogeneous in the brane direction x, for each type of 
perturbations, we can apply the Fourier transformation with respect to x — (x p ) to the 
perturbation variable as 

Sg liv = h lll/ {t,r i z i )e ik ' x . (5-4) 

Since the background metric is static, we can further apply the Fourier transformation with 
respect to t to h^ u if necessary and assume that 

V « e ' iult - ( 5 ' 5 ) 

Hence, we can reduce the Einstein equations for perturbations to a set of ODEs with respect 
to r. In this section, we assume that Jff n is compact. 

5.2. Tensor perturbations 

The equation for tensor perturbations (13-341) with tt = reads for the present system 

- %H T + ld r (r n fd r H T ) - f ($±*IL + A: 2 ) H T = 0. (5-6) 
Let us define the energy integral for a tensor perturbation by 

POO 

E := / drr 
Then, from the perturbation equation, we have 

oo 

(5-8) 

If there exists an unstable solution Ht oc e~ luJt with Imw < 0, it must fall off exponentially 
at r —>■ oo and vanish at the horizon from the above equation, provided that the solution is 
uniformly bounded. For such a solution, E becomes constant and contradicts the assumed 
exponential growth because all terms in the energy integral is non-negative definite. Hence, 
the black brane solution is stable for tensor perturbations. 

5.3. Vector perturbations 

5.3.1. Basic perturbation equations 

Basic gauge-invariant variables for vector perturbations are given by F a (t,r) with a = 
t,r,p (p = 1, ••• ,m). Among these components, we decompose the part parallel to the 
brane, F p , into the longitudinal component Fk proportional to the wave vector k p and the 
transversal components F^ as 

F k = ik p F p = d p F p , (5-9a) 
F^ = F p + ^F k . (5-9b) 



1 T 



f{H' T f + 



2K 



k m 



(5-7) 



E 



r n fH T H' T 
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With this decomposition, the perturbation equations ( 13-48aj) and f l3-48al) can be written 
as the four wave equations 

iflf F, - + " / + 7 2 + A V , = (/' - ^) ftF„ (540a) 

j3, 2 F' - -L-^ir"-^) + 2("-l)/+m„+<;V 2 F , = / ^ (J 1Qb) 

W - IflKrVflW + r/ ' + " /+ 2 m " + feV F t = -F', (540c) 

- f d 2 (F^/r) - -^d r [r" +2 /3,(F»] + (k 2 + ^f) (F x /r) = 0, (5-10d) 
and the constraint 

- J d * F t + ^=i d r( rn ~ l f F r) +F k = 0. (5-11) 

With the help of this constraint, the second of the above can be also written as 

- JLa r (r-MF) + ("-DPz-rn + m. + f^ = } , Fi (542) 

Clearly, the transversal part Fp 1 decouple from the other modes and each component 
obeys the same single wave equation. Further, each of (Ft, F r ) and (F r , Fk) obeys a closed set 
of equations, and the remaining components Fk and F t , respectively, are directly determined 
from them with the help of the above constraint equation. 

5.3.2. Master equation 

Let us take F r and Fk as fundamental variables and set 

W-=( ) (5-13) 

\(n + l)r n / 2 ~ l F r + r n / 2 F k J ' 

Then, the perturbation equations can be put into the form 

= (-D 2 + V + fA) V, (5-14) 

where V is the scalar potential 

V = f 



m v ^ 2 _|_ n ( n + 2) 

ry 2 /j^yj 2 



(5-15) 



and A is the matrix potential 

/ 2fc 2 _,_ (n+2)f 2fc 2 
_ [ n+1 ' 2r n+: 

n r 

*?.r J 



A =\ n+1 ^ ^ * „ I • (5-16) 



n+1 n+1 It • 
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In order to see whether this set of equations can be reduced into decoupled single equa- 
tions, we introduce a new vector variable <P by 



cp = Q& + p&\ 



(5-17) 



where P and Q are matrix functions of r that are independent of uj. If we require that <5 
obeys the equation of the form 



<P" + (to 2 -V - W)<P = 



(5-18) 



with a diagonal matrix W, we obtain constraints on V and B. 

For the exceptional mode with m v = 0, these constraints are satisfied, and we find that 
for the choice P = 1 and 



Q 



k^r n+2 r k^r 

m-4-1 9.r J n+1 

k 2 r I n £ 
n+1 ' 2rJ , 

W is given by the diagonal matrix whose entries are 



n+1 2r 
k 2 r 
n+1 



The corresponding equations for decouple to 



&! + (to 2 - Vi)$i = 0, 
V l = f 



k 2 + n + 2 



v 2 = f 



4r" 



n 



2(3n + 2)M\ 
n + 4: i — —1 — 1 



„n—l 



2nM\ 



V 2 is clearly positive. Further, in terms of the S-deformation with 



V\ is transformed into 



2r J 



V 1 = k 2 f > 



Hence, this system is stable for this exceptional mode. 

If we apply the same transformation in the case m v ^ 0, we obtain 



2m v fh 



r(r 2 oj 2 — m v f) 



fd r + co 2 -V 



(p 



fh 



(n + l)(r 2 cj 2 — m v f) 



B<P, 



(5-19) 



(5-20) 



(5-21) 
(5-22) 

(5-23) 



(5-24) 



(5-25) 



(5-26) 
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where 



(n + l)M 

h = l ~ { —^-' (5-27) 



V = f 



m v 2 n 2 + 2n + 4 (n 2 + An + 2)M m v fh 



+ k Z + — : - + 



r 2 4r 2 2r n+1 r 2 {r 2 ui 2 — m v f) 



(5-28) 

. (n + l) 2 oo 2 + 2m v k 2 -2m v k 2 \ 

B = \ y ' 5-29 
1 2m v k 2 -{{n + lfuo 2 + 2m v k 2 }) V ' 

Since B is a constant matrix with eigenvalues 

A = ±{n + l)u [(n + l) 2 tu 2 + Am v k 2 ] 1/2 , (5-30) 

we can reduce the set of equations for <P to decoupled single second-order ODEs. However, 
these equations are not useful in the stability analysis because their coefficients depend on 
uj 2 nonlinearly and have singularities in general. 

5.3.3. Stability analysis 

Since we cannot find a convenient master equation, let us try to analyse the stability by 
directly looking into the structure of the set of equations (15-141) . The subtle point of this set 
of equations is that the operator on the right-hand side is not self-adjoint because A is not 
a hermitian matrix. Therefore, we cannot directly conclude that uj 2 is real. 

Allowing for the possible existence of the imaginary part of u 2 , we obtain the following 
two integral relations from the above equation: 

f°° dr 

Re(uo 2 )(V,V)= / -[(W 1 ) 2 + ( J D^ 2 ) 2 + /f/ 1 |^ 1 | 2 + /[/ 2 |^ 2 | 2 ], (5-31a) 

Jr h J 

4k 2 f°° 

lm(u 2 )(&,&) = / drlm^&z). (5-31b) 

Here, D = fd/dr and 

Ui = ™v + ^+l k 2 + n(n + 2) + (n + 2)f 
r 2 n + 1 4r 2 2r 

4k 2 n + 1 f 

U 2 = U X --^-- ^-±L (5-32b) 
n + 1 r 

By applying the S-deformation with S = f-f to &2, the right-hand side of the equation 
corresponding to Re {uj 2 ) is deformed to 

771 71 — 1 

DV 2 ^{D + S)V 2 , U 2 ^^ + ——k 2 . (5-33) 

r A n+1 



*' In Ref 38 ' the author derived a well-behaved single master equation of 2nd-order for the black string 
background. There, the author took the gauge in which f z = and Ht — 0. Such a gauge cannot be 
realised in general because the gauge transformations of f z and Ht are given by 8f z — ~Sd z (L/S) and 
5Ht = k v L/S. If we set f z = 0, we cannot change the ^-dependence of Ht in general. 
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Therefore, if we assume that u 2 is real, as is assumed in most work, we can conclude that 
the system is stable against vector perturbations. However, we cannot exclude the possible 
existence of an unstable mode with Im (u 2 ) 7^ 0. 

5.4. Scalar perturbations 

5.4.1. Perturbation variables 

The gauge-invariant variable set F a b in the general formulation can be decomposed into 
the scalar, vector and tensor parts by their transformation behavior with respect to the 
brane coordinates as 

Scalar part: F tt , F tr , F rr , F kt , F kr , F kk , F±. 
Vector part: F ±pt , F ±pr , F ±pk . 
Tensor part: F± p ± q . 

Here, 



F ka = d p F pa /(ik) = (k p /k)F pa , (5-34a) 

F± pa = F pa - {k p /k 2 )k q F qa = F pa - (k p /k)F ka , (5-34b) 

F kk = {Wjk 2 )F w (5-34c) 

F ± = Fl - F kk , (5-34d) 

F_L P ±q = F Lpq - {k q /k)F ±pk - -^—F L (5 pq - k p k q /k 2 ). (5-34e) 



The remaining gauge- invariant variable F in the general formulation also belongs to the scalar 
part. Note that the vector and tensor parts do not exist for the black string background. 

5.4.2. S-mode 

First, we consider the exceptional mode with k s = 0, which is often called the S-mode. 
For this exceptional mode, the general gauge-invariant variables reduce to F^ = f a b and 
F = Hl due to the non-existence of corresponding harmonic vectors and tensors. These 
variables are not gauge invariant and subject to the gauge transformation law 

f 

5H L = --T r , 5-35a) 
r 

5f tt = 2iuT t + ff'T r , 5f tr = tujT r - f(T t /f)', 

5 f rr = -2T r ' - (f'/f)T r , (5-35b) 

5f tk = iuT k -ikT u 5fr k = -U-ikT r , 5f kk = -2ikT k , (5-35c) 

Sf± pt = iuT ±p , 5f Lpr = -T' ±p , 5f ±pk = -ikT ±p , (5-35d) 

Sf± = 0, 5f ±p±q = 0. (5-35e) 
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In particular, we have 



6 (ftk + ^fkk) = -iW, 5 (j rk + = -iW- (5-36) 

From these, we can construct the following five gauge invariants for the scalar part: 

r 2-n X = _ F± ^ (5-37a) 
r 2 ~ n (X -Y)= F r r - 2rF' - ^ - 2^ F, 



(5-37b) 



r 2-n z = pr + _±_ {F ty + _ ^_ ^ F ) ? (5 . 37c) 

,2-TlTT-t _ pf , ^ ( jpt r f' tA 10 



r v ~ 2tj i F ' " t F J " S7 F "' (5 ' 37d) 



r 2 -l/' = *J - ikrF + 'If^ (5-37e) 
For the vector part, we adopt the following two gauge invariants 

r 2 - n W l p = F{ p - ^F ±pk , (5-38a) 

r 2-n w r = ^ + (5 . 3 g b) 

(1) Tensor part. 

First, we study the stability in the tensor part. The perturbation variable of this part, 
F± p ± q , follows the closed equation 

- f(r n fF' ±p±q )' + (k 2 f - cv 2 )r n F ±p±q = 0. (5-39) 

From this we obtain the integral relation 

roo poo 

: / dnr n \F ±p±g \ 2 = / drr n [f\F' ±p± f + k 2 \F ±p±q \ 2 ] - [r n fF^F> ±p±q \™ . (5-40) 

Jr h Jr h 

If there exists an unstable mode with 

uj = uj 1 + iu 2 ; u>2 > 0. (5-41) 
a solution that is bounded at the horizon behaves as 

F ±p±q ~ e~ iu,r * (5-42) 
near the horizon. Next, at infinity, the solution behaves 

F ±p±q ~ - 1 -L^Z u (Vuj 2 - Wr) ~ r- n / 2 exp(±2V^ 2 -A; 2 r). (5-43) 



u 2 
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Therefore, for an unstable mode that is uniformly bounded, the boundary term in the above 
integral relation vanishes and the integral at the left-hand side converges. This implies that 
u 2 > and leads to contradiction. 

(2) Vector part. Next, for the vector part, we obtain the following two equations for the 
gauge- invariant variables W p and W p : 



-iuo 



n-2 , /' 
/ 



+■ - 1 w± p 



+ 



k 2 f - LU 2 



W r ±p = 0, 



-iuWi p + Wl p ' + -W r ±p = 0. 



(544a) 
(544b) 



Therefore, we can set 

Wl p = r~ 2 <P, iuW{ p = r~ 2 <P', (545) 
and the perturbation equations can be reduced to the following single master equation for 

- f(r' n f<P')' + (k 2 f - uo 2 )r- n $ = 0. (546) 

By the same argument for the tensor part, we can show that this equation does not have a 
uniformly bounded solution with Im (u) > 0. 

(3) Scalar part. 

Finally for the scalar part, the perturbation equations gives the closed lst-order set of 
equations for X, Y, Z, V*, 



X' = 



k 2 rHf 2 



r uj -co <i k r f + n — n(n + \)x + 

1 



2 + 



-J^(nui 2 -k 2 H)Z 



X + 



k 2 rfH 



nuj z 1 



3n 2 + 2n - 1 
I 
n+1 



-x 



x ) + k 2 H 2 \ Y 



Y' = 



k 2 f 2 H 
1 

k 2 rf 2 H 



krfH 



{2uj 2 r 2 + (n- l)xH) V 1 



(547a) 



r V - uo 2 { 2k 2 r 2 f + n - n{n + l)x + 



3n 2 + 2n - 1 



-x 



+r 2 k 4 f 2 - In - (n 2 + l)x + 



n + 1 



{n + If 



x 2 k 2 f 



X 



k 2 rfH 

2n - : ^-k 2 f)z + J 



noj I 1 x 1 + (n — l)k In x + 



5n 3(n + l) 



-x 



Y 



k 2 f 2 H 
. (n — l) 2 x 



krfH 



{2rV - 2k 2 r 2 f + (n - l)xH) V\ 



—i- 



2ujt 2 



X + {r 2 uj 2 + {n - l) 2 x}Y 
2r z u 



rf ' 1 



n + 1 



x ) Z + ikfV\ 



(547b) 



(547c) 
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(V*)' = 



CO 



r 2 u) A + uj 2 < 2k 2 r 2 f + n - n(n + l)x + 



3n 2 + 2n - 1 



-x 



2rk 3 f 2 H 

i 



k 3 f 3 H 
1 



, , n + 1 \ l2 J , 1N 3n 2 + 2n- 1 
- //a,' - ! 1 - —7 — ./' J + k in n[n + l)x H 



x 



+ 



r uj + oj \ k r j 



2 

n — 1 



rk 2 f 2 H y ' ~ ' ~ V" " 2 

and the expression for V r in terms of these quantities, 



xH)+(n- 2)k 2 f 2 H } V\ 



V r = 



i(uj 2 + k 2 f) 
2nk 3 rf 2 



{(uj 2 + k 2 f)r 2 + nf 2 } X + nf 2 Y + 



2inuj r 
uj 2 + k 2 f 



Z + 2ukr 2 fV t 



Here, x = 2M/r n ' 1 . 

From these equations, we find that X obeys the closed 2nd-order ODE 



-f(fXy+(n-A)£-X' 
r 

-u; 2 + f(k 2 + ^{l + (n-2)x} 



+ 

which can be put into the canonical form in terms of <P defined by 

X = r n/2 ~ 2 <P 



X = 0, 



as 



f(f&)' + [-uj 2 + / {A; 2 + £^(n - 2 + nx)}] <P = 0. 



It is clear that this equation does not have an unstable mode. 
Next, let us define the new variable Q by 

Q := PX + nf ( 1 - Y + 2mu;rZ + 2kcur 2 fV t ; 



P :-. 



n + 1 (n — l)x ( n+1 
■x _ , „ „ n — x 



2n ' 2k 2 r 2 \' 2 

n + n(n + l)x - (3n 2 + 2n - l)x 2 , 



o o 1 2 2 

u; r H — - — r 



2n 



Then, we find that Q satisfies a closed 2nd-order ODE mod X = 0: 

- f(fny + A] a + (-^ 2 + v n )n = bx, 
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where 



k=0.5_ 





Fig. 8. The effective potential for S-modes 



A 



f 



[{An 2 + 2(n + l)(n - 2)x) k 



2 2 

r 



Ar 3 gH 

+n 2 (n - l)x {3(n + l)x - 2(n + 2)}] 



(5-55a) 



V, 



f 



n 



[2 {2n - (n + l)x} 2 k 4 r 4 



8r 4 g 2 H 

+ {8n 2 (n + 2) - 4n(n + 2)(3n 2 + n + 2); 
+2n(n + l)(8n 2 + 5n + 5)x 2 - (n + 2)(3n - l)(n + 1)V} £; 2 r 2 
+n 2 (n - l)ac + l) 2 x 3 - 3(3n - l)(n + l)a; 2 
+4(2n 2 + 2n - l)x - An 2 }] , 



B 



fg 

nr 2 H 



{(n + 1)uj 2 - k 2 } [-2£; 2 r 2 (l - nx) + n{n - l)x(n - x)] 



^ - * + 2r 2 

By the transformation 



x, g := n — 



n + 1 



-x. 



Q = r n/2 gH&, 
we can put this equation into the canonical form 

- f(fW)' + ( — uj 2 + = r n/2 gHBX, 

where 



(5-55b) 
(5-55c) 

(5-55d) 
(5-56) 

(5-57) 
(5-58) 
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U = k & + ( n + 4 J k \ n + 2-3nx) 

" ( " {3n(n + 2)-(2n 2 + 3n + A)x}x 



4r 4 

n 3 (n — l) 2 9/ „ . . . 

+ ^ -^x 2 (n - 2 + nx). (5-59) 

16r b 



This potential has a deep negative region for < k < k n with some constant k n dependent 
on n. It has been shown by numerical calculations 37 -*' 38 - ) that the eigenvalue u 2 becomes 
negative for some range < k < k c , as first pointed out by Gregory and Laflamme using a 
different reduction. 32 ) 

5.4.3. Generic scalar perturbation 

(1) Tensor part. The tensor part of generic scalar-type perturbations obeys the decoupled 
2nd-order ODE 



" f(r n fF' ±p±q Y + (-u 2 + k 2 f + r"F ±p±q = 0. 



(5-60) 



It is obvious that this equation has no unstable mode. 

(2) Vector part. 

In terms of the gauge-invariant fundamental variables 

V t = r n ~ 2 F t±p , V r = r n ~ 2 F r±p , V k = r n ~ 2 F k±p , (5-61) 
the perturbation equations for the vector part are expressed as 

-(fV r )' - iuf^Vt - ikV k = 0, (5-62a) 



.00 



2 + V r + ik (v k ' - = 0, (5-62b) 



-r-V-VW + (-y + ^ + n ~^f + n ~^f) v k 

~V t + ik f(fV r y + ^-V r ) = 0, (5-62c) 



V r 



{^T 1 + Vi ^r) v * ~ ( V r + l V r) + >«*V k = 0. (5-62d) 
By eliminating V t and introducing the new variables Y and Z by 

<p = f Z y J ; ^ = -^W 2 - 2 z, v; = /-V"/ 2 - 1 ^ (5-63) 
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this set of equations are reduced to a set of two ODEs, 



D Z <P -{-u 2 + k 2 f + '—^f + 



n + m n(n — 2) 



f}<P = A<P; 



-2kf 

v J It ^ t 



2r 



(5-64) 
(5-65) 



This set of equations has the same structure as that for vector perturbations and can be 
shown to have no unstable mode if uj 2 is real. 

(3) Scalar part. 

Finally, we discuss the scalar part of the generic scalar-type perturbation. Utilising one 
of the Einstein equations 

E T = 2{n - 2)F + Fl = 0, (5-66) 
the basic perturbation variables can be expressed in terms of X, Y, Z,V t ,V r ,S and & as 



F\ = X + 2F - k 2 fV\ F r r =Y + 2F, F[ = ituZ, 

F r k = ikV r , F* = ukV\ F kk = S + uj 2 V 1 + 2F, 

2(n + l)F = -V - X - Y - S - (uj 2 - k 2 f)V\ F ± = V. 



(5-67a) 
(5-67b) 
(5-67c) 



Here, Q = r n ~ 2 Q in general. 

In terms of these variables, the Einstein equations can be reduced to the decoupled single 
equation for 



-co 2 + k 2 + 



n + m 



f 



- r- n f(r n f&')' + 
and the regular lst-order set of ODEs for X, Y, Z, V\ V r and S, 
Z' = X, 



W = 0. 



X' = ^X + 



+k 2 f'V t , 



f 2\ r 1 / lu 2 ;2 m + n\ „ 
^Y + -{ — - + k 2 + — 1 Z 



f rj f\ f 



Y' = f f (X-Y) + f 2 Z+j(v r - f -fv t ), 

(v r y = -s, 

S' - T —S--Y-„:-\ 



lf\rt , 1 ,2 n + m 



f f\f 



— k — 



(5-68) 

(5-69a) 

(5-69b) 
(5-69c) 
(5-69d) 
(5-69e) 



k 2 r 2 ff 2 (Vy 



2u r + (n — l)x I n — 



n+1 



-x 



X 



+ 



2u 2 r 2 - 2{k 2 r 2 + n + m)f + 2n - 4nx + ^ + ^ ~ 2 



-x 



Y 



39 



+ - [— 2nuj 2 r 2 + (n — l)x(k 2 r 2 + n + m)] Z 



(n- l)k 2 x 2 



n — 5 



a: )fV t - 2k 2 r n 



n + 1 



x V r 



-2k 2 r 2 fS, 

where x = 2M/r ri ~ 1 . 

If we define X 1: X 2 , X 3 by 



(5-69f) 



X X = Z, 

X 2 = r(X + Y)-nZ + k 2 rfV\ 
rf 



X 3 = - 1 + 



2nf 



1,2 2 

r(X + Y)- nZ] - ^-f'V\ 



2n 



and introduce $ by 



/ r -n/2 Xi \ 



(5-70a) 
(5-70b) 

(5-70c) 



(5-71) 



r -n/2+l / -l X2 

we can reduce the above set of lst-order ODEs to the set of 2nd-order ODEs of the normal 
eigenvalue type as 

uj 2 <P = (-D 2 + V + W)$. 

Here, V is the scalar potential 

/ 



(5-72) 



4 r 2 



[4(m + k 2 r 2 ) + n 2 - 2n + n(n + 4)x] , 



and W is the following matrix of rank 3: 

(n 2 - l)xf 2 



W n = 0, Wi 



12 



2f 



13 



w 2 i 

W^22 
W 3 1 



{4 - 2{n + l)x} k 2 r 2 - 2(n - l)mx - n(n 2 - l)x 



rf 



nf 2 

- J 2 ~, W 23 = 0, 



1 



2nr 2 f 



[2(n - l)x(n - 2 + x)r 2 k 2 + {An + 2n(n - 5)x + 2(n + l)x 2 } m 



+n(n + l)x [2n 2 - (2n 2 + 3n - l)x + n(n + l)x 2 }] , 



W, 



.32 



{n 2 - l)x{2 - (n + l)x}f 
2nr 3 
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(n + l){2-(n + l)x}f 



(5-73) 

(5-74a) 
(5- 74b) 

(5-74c) 

(5-74d) 
(5-74e) 



The original variables X, Y, Z, Vt, Vr, S can be expressed in terms of Xi and DXi as 

X = X[, Y = -X[ + -X 1 - { -^^X 2 --X 3 , Z = X X , (5-75a) 

r 2nrj r 
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(n + l)x-2n)„ , 1 



^ = - W (5 ' 75b) 



(n — l)x . / . 2r 2 k 2 + 2m + n(n + l)x 
2nk 2 r 2 ~Wr 3 2k 2 r 2 ' 



(„- 1), {(„+!),- 2}^ _ / A ._ 



2nk 2 r 2 f k 2 r 2 
2k 2 r 2 + 2m + n(ra + l)x , n/ . (n — l)x(2k 2 r 2 + 2m + n + n 2 ) 
= 2Pr 2 1 + kh^ 3 + 2fc 2 r 3 / 1 



-cu 2 r 2 + f(n + m + fc 2 r 2 ) 
Pr 3 / 



{n - 1)x x 2 + x 3 



2nf 



(5-75d) 



Hence, the equation for $ is equivalent to the original 1-st order system. 

Unfortunately, it is not possible to analyse the stability of this 2nd-order system by an 
analytic method, partly because it is not a self-adjoint system. However, all numerical cal- 
culations done by various authors have found no evidence of instability for this system. 37 )' 38 ) 

§6. Summary and Discussion 

In this lecture, we have explained the gauge-invariant formulation for perturbations of 
a class of background solutions to the Einstein equations that include various practically 
important spacetimes as special cases. Then, we have illustrated its power by applying it to 
the stability problem of static black holes in higher dimensions and flat black branes. 

These two systems have one important common feature in addition to staticity that 
the background spacetime is of the cohomogeneity one. That is, the isotropy group of 
the spacetime has orbits with codimension one, and roughly speaking, the spacetime is 
inhomogeneous only in one direction, say r. In this case, the perturbation equations for the 
system can be automatically reduced to a set of ODEs for functions of r with the help of 
the harmonic expansion. This applies to a rotating black hole case as well 56 '' 

There exists however one crucial difference between the two systems. In the static black 
hole case, the perturbation equations can be reduced to decoupled single 2nd-order ODE, 
and the stability problem is formulated as an eigenvalue problem for the corresponding self- 
adjoint operator. In contrast, in the black brane case, it appears to be impossible to reduce 
all the perturbation equations to decoupled single master equations. Further, the eigenvalue 
problem for the stability issue cannot be put in the self-adjoint form even if we allow for a 
mult i- component expression, except for some special modes. Nevertheless, numerical calcu- 
lations indicate that there exists no eigen-mode with an imaginary frequency 37 -*' 38 - ) There 
must be a profound reason behind this result. 

The gauge-invariant formulation developed in $3] can be also applied to spacetimes whose 
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cohomogeneity dimension is greater than one. Although this implies that the formulation can 
be applied to perturbations of rotating black holes in higher dimensions such as the Myers- 
Perry solution 57 ^ and its generalisation to non-vanishing cosmological constant, 58 '' it may 
not be practically useful in most case, because we obtain a couple set of partial differential 
equations in a reduced spacetime with dimensions smaller than the original one. However, 
in some special cases, we obtain a single PDE that is separable to ODEs. For example, for a 
Kerr(-adS) black hole that rotates in a two-dimensional plane, we can classify perturbations 
into tensor, vector and scalar types as in the static case, and among these, the perturbation 
equation for the tensor-type perturbation is separable. 43 ** There may exist other cases in 
which similar phenomena happen. 
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